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Abstract 

A generalization of the addition relation for the Riemann theta functions and its limiting version for 
exponential functions appearing in soliton type equations are reported. The presented form seems to be 
particularly useful when processes in A'^ + 1, (TV > 1) space-time are analyzed. The commonly applied 
bilinear and trilinear approaches, restricted to the pure soliton processes, represent particular cases of the 
reported formalism. As an example, the dispersion equation for either quasiperiodic or soliton processes 
following 2+1 Calogero-Bogoyavlenskij-Schiff equation is derived. 

PACS: 05.45.Yv, 02.30.Jr, 02.30.Gp. 

Running title: Generalization of addition property... 



1 INTRODUCTION 

As it is well known, solitons represent the localized excitations and appear in the numerous branches of 
physics when nonlinear description is essential, as in superconductivity, plasma physics, fiber optics, domain 
wall or fluid dynamics, protein chains etc. Unfortunately, the well developed soliton theory practically deals 
with the one-dimensional dynamical systems. The same conclusion relates to the quasi-periodic processes. 
Because of the physical reasons, many research groups make efforts in order to find the methods for mul- 
tidimensional solitons analysis. The results obtained hitherto have however only a contributory character. 
Therefore below we present, it seems, a slightly more general approach to the truly multidimensional solitons 
and quasi periodic-processes. 

In order to illustrate the problem we start with a simple example of the famous Kadomtsev - Petviashvili 
equation 

iuyy = [4ui - &uux + u^xx]^ , (1) 

repeating some arguments of Dubrovin j^]. Symbols (.)^ or Ux denote relevant partial derivatives. Looking 
for the solution in form 

M = -u (x,?/, i) = -21nr (zi, ...Zg) , Zi = hx + ky + Wit + ZQ^i, i^l,...,g, (2) 
a substitution to the equation leads to 



^'^xxx^x ^~ 3 (^Txx^ ^ ^^xt^ '^x'^t ^^yy^ ('^y) 

if the first constant of integration is zero, and the second one is equal to — Sd. According to the Hirota 
approach, there is introduced the bilinear differential operator Dj. which if applied to the ordered pair of 
functions / {x, y, t) and g [x, y, t) is defined as Q 

if o g) := {d^f [f {x + w, y, t)g{x- w, y, t)]^^, . (4) 
Then equation can be written as 

+ ADxDt - 3 {Dyf \ (r o r) + Sdr^ = 0. (5) 
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Next, for d = 0, making use the standard Hirota procedure one can find multisoliton solutions of (|l|). 

It is worth notice that the case d ^ is not considered in frame of the bihnear technique since it leads 
to quasiperiodic solutions and periods of the periodic subprocesses tend to infinity as d — > 0. 

There is an opinion that a huge success of the bilinear formalism in soliton theory Hj' |^ can be linked 
with the addition property for r-functions. 

t(z + w)t(z- w) =^W^(w;e)Z(z;£), (6) 

where z, w G C^, r : ^ C, W, Z : x Z,^ ^ C , the sum is over a finite set and Z (z; e) enumerated by 
e GZj form a set of linearly independent functions. Exponential functions appearing in multisoliton solutions 
have this property as well as the Riemann theta functions leading in turn to the quasiperiodic solutions. 

The essential feature of the equation is that the product of shifted r-functions admits factorization. 

Relation (^) should be considered only as the t -function property. This means that putting either w = 
or z = we have (z) — W (0; s) Z (z; e) or r {w) r {—w) = J2e ^ ^) ^ (^i ^) ' respectively. If the 
proper determinants do not vanish, functions W (w; e) and Z (z; e) can be expressed by r-functions. 

(Nota bene, if r-function is identified with the Riemann theta functions 9 (z) - see below, for £,e'eZ2 , 
we have 9 {z + e' /2) = 9 {z - e' /2), then W {e'/2; e) = Se,e' - Kronecker symbol and Z (z; e) = 6*2 (z + e/2) .) 

Thus it seems that relation m) is really very important and there are a few practical reasons for such 
statement which we shall discuss below. 

First of all, without any additional assumptions for function having addition property one can easy 
calculate even derivatives of the equation (^, obtaining some combinations of derivatives of r— function 
logarithm, Q 

e 

(7) 

Ujki + 2 (3 X L,jLki) = Wmw,w,w, (0; e) [Z (z; e) /r^ (z)] , 

6 

(8) 

(0;e) [Z(z;e)/r2(z)], 

e 

(9) 

etc., where L := Inr (z), Li :— S^. Inr (z) and we use a shorthand notation (3 x LijLki) LijLki + Ln^Lji + 
LiiLkj i.e. including all permutations. 

Note that the "basis" functions [Z {z;e) /t^ (z)] for all operators are the same. Therefore, the differenti- 
ation rules (^, (H), (H), etc., automatically reconstruct the Korteweg - de Vries - Kotera - Sawada hierarchy 
giving a tool for a derivation of relevant dispersion equations [Q ' Q . 

In order to elucidate problem of dispersion equations, let us observe that (^) can be written as 

Lxxxx + 6 (Lxx) + 4Lj;4 — -^Lyy — ^8d, (10) 
and, due to the linear dependence of arguments of L = Inr (zi, ...Zg) on space and time coordinates, as 



Y kpkqkrh [Lz^z^z^z, + 2 (3 X L z^z,L z,.z,)] + ^ {'^kj.Wq - Slplg) Lz^z, = Sd. (11) 
p,q,r.s—l p.q—1 

The shorthand notation 3 x ... was explained before. Now, if r— function has the addition property, we can 
apply rules (|^) and (||) obtaining 



E 



q r s—1 ^p^Q^T-ks^^^WpWqWrWs (S^: ^) 

+ EL=i '{"^kpiu, - 3lplq) W^^^^ (0; e) + 8dW (0; e) 



Z{z;e) = 0. 



(12) 
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since (z) — W (0; e) Z (z; e) . Finally, the Kadomtsev - Petviashvili equation (|l|) has a solution if in 
a relevant class of r— functions (exponential or Riemann theta function), for any e ^1^% system of 2^ 
algebraic dispersion equations 
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(0; e)+Y, (4^p^^9 - ^Vi) (0; e) + MW (0; e) = (13) 

p,g,r,s— 1 1 

has a nontrivial solution. The broader discussion, particularly concerning quasiperiodic solutions, one can 
find in [Q^ ||^' . Here, we want to underline only a formal similarity of (^ and ([l3| ) equations and to 
underline that equation ( |l3|) remains valid also for d ^ 0, i.e. also for quasiperiodic solutions 

The second argument relates to another class of equations which are solved by means of so-called trilinear 
operator. The trilinear operator T and the complex conjugate trilinear operator T* are defined as 0' 

T{fogoh)^{d.,+jd,,+pd,,) /(zi)5(z2)/l(z3)Ui=..=.3=-, (14) 
T*{fogoh)^{d,,+fd,,+jd,,) /(2i)5(^2)/l(^3)Ui=..=.3=., (15) 

where j — exp (i27r/3) . 

Its application we illustrate by example of the Satsuma equation ^ 

^xx^yy-^ -^xx {-^y) ^yy (-^a;) ^^ '^^xy^x-^y (-^a;y) — 0; (1^) 

which then can be written then as 

{TxT^TyT; - T,X^) {FoFoF) = Q. (17) 

Observe now that the addition property for _F— functions in version is not applicable, since ( p^ 
consists triads of F derivatives in contrast to (||) where we have dealt with the pairs only. Therefore the 
left hand side of a relevant version of addition property (||), if exists, ought to contain the product of three 
shifted functions, (or even more). 

The next argument relates to discrete equations and their reduction to dispersion equations. For example, 
considering the KdV- type completely integrable difference-difference equation 



u {x,t + d) — u {x,t — d) = u{x,t + d) u{x,t — d) [u{x + d,t) — u{x — d,t)] , (18) 
u{x,t) -.^ t{x + d,t)T {x - d,t) /T{x,t + d)T {x,t - d) ~ 1, (19) 

it is seen that equation (^ enables to write the products t {x + d,t) t {x — d, t) and 

T {x,t + d)T (x, t — d) in a, compact form leading to the dispersion equation The shift w = d plays a role 
of a step in a difference equation. On the other hand, if the multidimensional version of difference-difference 
equation is considered, it requires an introduction of a few and independent steps with respect to each 
coordinate. This in a natural way leads to the form of addition property, where on the left had side the 
product of a few shifted functions appear. 

Now, one can ask about the class of functions having the property (|6|). The oldest and the most known 
functions having this property are the Riemann theta functions ' [Q ' ||] . 

e iz\B) = J2 exp (2 (z, n) + (n, Bn))] (20) 

where z G C^, B G C^^^ is the Riemann matrix, (i.e. symmetric with positively defined imaginary part), 
(z, n) := ^j^j- The equivalent of (||) takes then the form 

6i(z + w|B)6i(z-w|B) = (21) 
[exp [iTT (2 (z, e) + (e. Be))] 6* (z + Be\2B) exp (i27r < 2w, e >) 61 (w + Be\2B)] 
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where symbol e G denote g-fold sum over ej=0,l;i — l,...,g. This relation can be also written down 
by so called - functions with characteristics. 

The second class is represented by exponential functions 



E z\B 



= J2 



in ( 2 (z, n) 



n, En 



(22) 



which appear in solutions of standard soliton equations. Matrix B E C^^^, (although sometimes it is 
convenient to assume that diagonal elements of B are real). In our opinion this is just a source of the 
bilinear operator applicability. The combinations of theta and exponential functions leading to processes 
when solitons propagate on a quasiperiodic background form the next class of functions having addition 
property |]. 

The fourth class, which to our knowledge has no physical application, is generated by integral of product 
of Gaussian function and Riemann theta functions with characteristics, where integration takes place with 
respect to the second characteristics. 

Concluding our motivation, the interpretation of the shift w in (|^) as a step suggests that for solitons in 
+ 1 space-time a few independent steps Wi, wn should be introduced. This leads to the generalization 
of the relation ( ^l| ) or more precisely of equation (^ and such generalization is just the aim of this note. We 
believe that the reported below relations can be useful for some multidimensional soliton type problems. 



2 THE GENERALIZED ADDITION PROPERTY 

In order to repeat the procedure as for (|^-|^) in case of a few independent variables it would be sufficient to 
have relation 



(z + u("))t(z + uW) ...t (z + - 



.(1) 



W 



e Z(z;£) 



(23) 



eC9, j = 0,...,J-l. 

,(j) 



for some class of r -functions, where u'^' — u'^^^ (w^ 

The reason is following: the derivatives of ( p^ ) with respect to Wj"" will relate only to 11^— function, 

leaving Z (z; s) unchanged. On the other hand the derivatives of l.h.s.of (p3| ) with respect to wj"''' one can 
change into derivatives of r— functions with respect Zi which is necessary for an application in the soliton 



theory since usually z is linear in space and time variables, (z := li^x + 



cut e CS). Since 



as a r -function, either exponential or Riemann theta functions are usually chosen and since exponential 
functions can be considered as a particular case of the second ones, we shall start from the Riemann theta 
functions. 

There are numerous transformations for 9 - functions, but to our knowledge the J-th order 

addition relation, as here, was considered only by Koizumi and cited in But his factorization of 
the r.h.s. was completely different than required in ( |23|) and thus it is rather useless for our purposes. 
Nevertheless for the standard 6 - functions according to ( |20[ ) one can prove the identity which coincides with 
our demand (|2|), 

Theorem 1. 

j-i 

If z,u('^) e Cf, fc = 0,.., J - 1, such that = 0, then 

k=0 



,(2) 



exp [iTT (2 (z,£) 



{e,Be))] 0{Jz 
I 



\b)...0(z- 

Be\JB) X 



u^'-'Mb] = 



(24) 



X exp 



< u 



(0) 



e > 



u(o) - u(i) 

U(0) - u(2) 



V u(0) - u('^-i) 



Be 
Be 

+ Be 



2B 


B 


B 




B 


2B 


B 








2B .. 




B 


B 


2B 


) 
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where all 6— functions are of order 5, with exception of the last one on the r.h.s. which is of order (J — 1) g. 
The proof we report in Appendix 1. 

Observe that the g-fold sum now is over Sj = 0, 1, J — 1, (it contains elements). If instead of J 
vectorial parameters u'-'^\ among which only J — 1 are independent, the new J — 1 vectorial parameters 
w^'^^ (fc=l,...,J — 1) are introduced by 



u 



(k) _ 



E'nJj-k+i j™w(™+''^-^) + j™w("+'=~i) , for = 1, J - 1, 



for k — 0, 



J-i 



(25) 



where j — exp (z27r/ J) , the relation ^^u^'^^ = is satisfied automatically since X]fe=i = ^ ^k=i 



k=0 



Thus equation (gj) with parameters given by (|5|) is satisfied for arbitrary set w^'^^ e C^. The quantity 
j is of course J-th root from unity. 

For fixed J, equations (E5h can be inverted, leading to 



rik) 



J 



= 1,...,J-1, 



(26) 



where 6k,j-i is the standard Kronecker symbol. Note that the choice of w parameters as follows from 



( p5| ) or (26) is not unique. This one adopted here however, gives a correspondence with trilinear operators 
introduced earlier in soliton theory, 

Applying a procedure denoted as the soliton limit [Q to the identity ( p^ ) we are able present rewrite 
relations (|2^) for exponential functions (p^). 

Theorem 2. 



= ^ |exp i7r2{^(z + u(°)^ ,£^ + (|£,i?e^ 
c (e, m) exp 

E 



«7r 2 ( ,/z + eS, m ) + J ( m, _Bm 



',7-1 



X > < C ( J, £, n) exp 

n(i),..,n(-'-i)eZ? I 



,7-1 ,7-1 



vfe=i 



fc=i i=k 



where "cut-off" functions c (e, m) and C ( J, e, n) are given by 

g 

c(£,m) = J]^ (J,„^.,o +(5™^,i(5e,,o) ; 



C(J,e,n) = n 



■j-i 



7-1 



( J) = E 



.fc=0 



E '^^3.'= i^'^AJ)'k-i (1 - 4,0) + SNj(j),k) 

sW=u(")-uW, /fc=l,..,J-l, 



fc=i 



(27) 



(28) 
(29) 
(30) 



u^*^) are defined by ( p5| ) and i? matrix is such that Diag Im B = 0, which however does not infringe a 
generality. 

The proof follows from the observation that exponential functions E (^1-^) can be obtained from theta 
functions 6 {z\B) as a limiting relation 
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e(z\b)= lim 9 ( z--iDeQ\B] , D := Diaglm{B) , B:=B-iD, (1, 1, 1)* , 

V / all Dii^oo V 2 / 



which was cahed - the sohton hmit, Q. 

Similarly as in the case of standard addition property one can calculate some combinations of 
derivatives of M^-function. They can be reduced always to the different combinations of derivatives of r 
-function logarithms, similarly as in case of the standard addition relation (Q-^. The sketch of procedure 
we present in Appendix 2. 

3 PARTICULAR CASES 

For J = 2, since j = —1, the identity ( p^ ) reduces to the commonly know n (pT| ) formula, and similarly for 
exponential functions. For J = 3, since then j = exp (i27r/3) , the identity (^4[) reduces to 



0(z + wW+w(')|Bj 0[z+jw(i)+/w(2)|Bj ^?[z+jVi)+jw(2)|Sj = (31) 
= W^(wW,w(2);£)z(z;£), 



M^(w(i),w(2);£j = (32) 
= exp fi27r < ( w^^^w'^M , e > 



(l-j)w(i) + (l-j2) w(2)+Be 
(l-j2)w(i) + (l-j)w(2)+Be 



2B B 
B 2B 



Z(z;e) = Y exp[in{2{z,e) + {e,Be))] ei3z + Be\3B) . (33) 

where z, w^^^ w*^^-* e and are arbitrary. This relations rewritten for exponential functions shows a close 
relation with trilinear operator analyzed in papers of . 

For J = 4, J = i and we obtain the product of four 9 - functions in version, which differs however from 
that reported by or |]. 

4 AN EXAMPLE OF APPLICATION 

Let us consider the Calogero-Bogoyavlenskij-Schiff (CBS) equation (r/ = 0) , which follows also from a reduc- 
tion of self-dual Yang-Mills equation and its modification by Grammaticos-Ramani-Hietarinta {rj = 1) , Q], 
of which soliton solutions were reported in Q 

+ ^'^xxxy + ^x'^xy + ^'^xx'^y + V^d'^^yyy = 0. (34) 

After substitution $ = 2Lx , where L = In t this equation reduces to 

V-^yyy ~^ ^^xxt ~^ ^xxxxy ^~ ^-^xx-^xxy ^" ^-^xy-^xxx (^^) 

We assume that the integration constant C is independent of y and t as well. Moreover, if t = r (z) = 
T (zi, Zq) , Zk = K.kx + VkV + i^kt, fc = 1 , ...,g , thc equation (|5|) becomes {Lijt := dz^dz^dz^ Inr (z) ) 



{fl^iVjVk + 4KjKjWfc) Ly-fe -I- (36) 

"pr X/ K^'^il^j^kKl'^m + 4:Kil^jKkKlKm) {Lijklm + ^^LijLklm)] = C 



ijk 

+ 12 



ijklm 
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As the T-function we can choose either 9 -functions or exponential E'— functions, (or even some combi- 
nations of both in spirit of those reported in Q]). 

On the other hand for J = 3, substituting w = w'^^^v = w^^^ in ( |3T| ) for typographic reasons and 
denoting 

S* (z, w, v) = T (z + w + v) r (z + jw + j^v) t {2, + j^w + jv) , (37) 
simple but a little tedious calculations give the necessary relations 

dwiWjw^S (z,w,v) |w=v=o/5' (0,0,0) = dy^y^-UkS (z,w,v) |w=v=o/'5' (0,0,0) = 'SLtjk (38) 
S (z, W, v) I 

w— V— 

/S (0, 0, 0) = 3 [L.jkim + l2U,Lkim] ■ (39) 

Derivatives of ( |3l| ) with respect to w and v deal with the W function and since S (z, 0, 0) = [9 {z\B)f' , 
we have (see Appendix 2) 

^*jfc = \Y. [^-.-.-'c lw=v=o Z (z; £)] (40) 



EG 

1 



Lijklm + ^'^LijLklm — 2 



VKu,,^,^u,^u,,t,„ (w,v;e) |w=v=o ■Z'(z;£) 



^3 



(41) 



Coefficients W(w,v;£:) and functions Z(z;e) are given now either by ( |32| ) or (|33|), respectively, if 
quasiperiodic solutionsare considered or by their soliton limit in spirit of (p7|). 

Considering functions ^ (z; e) as independent, (|6|) reduces to the system of 3^ algebraic dispersion 
equations for any e G Z3 

{'qViVjVk + 4KiKjUJk) Ww.WjWt (w, v; e) |w=v=o + (42) 

ijk 

KiKjK/cKiI^mVF^.-u)jmfc-u),i;„ (w,V;e) |w=v=0 = 3CVF (w, V; e) |w=v=0 

which determines relations between Ki^Vi^oji {i = I...5) and also C, if ry = 1. Its nontrivial solution, if 
exists, gives the quasiperiodic solutions of 2-1-1 CBS equations in form of $ = 2 [ln^^ (z|i?)]^ with — 
KiX + UiU + ujit. For higher genus it gives also the additional conditions on matrix B, since system ( ^2| ) 
can be overdetermined. then). 

Multisoliton solutions of this equation, reported in paper ||^ , were found by means of the trilinear operator 
formalism, where the third roots appear as above, but they can be also derived using equation ( [2^ ) taken 
for J = 3. There is an interesting similarity between trilinear soliton version of CBS equation reported in 
Q and the left hand side of equation (^) which when C 7^ is valid also for quasi-periodic processes. It 
follows from the fact that trilinear operators T and T* used in papers Q and coincide with used here 
operators d/dw and d/dv, respectively. 

In conclusion, it seems that the presented here generalized addition formula can be useful for other + 1 
soliton type equations. 
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6 APPENDIX 1 



The proof of Theoreml. 

Lemat 1 

2B .. B 
2B .. 
B .. 2B 

B) is also Riemannian. Indeed, both matrices B and B'aie symmetric. Moreover, if (n, Bn) > 0, for any 
n G (with equahty for n = 0), then 



If B is Riemannian then also matrix B' = 



( with diagonal blocks 2B and otherwise 



j-i 



'j-i 



(ni,...,n^_0,B'(ni,...,n^_0 ) = ^ (n^, Bn,) + ( ^nJ,B^nJ )>0, 



fc=i 



(with equality for all nk= 0) 



By the definition of 0- function (20), the l.h.s of (M) takes form 



J2 - J2 



H ... 

ri(0)eZs n(-'-i)eZ9 



.7-1 



,7-1 



[2 ( (z + u(0)) , n(") ) + (n(") , Bn(") )] + 
+2E/r/ ((z + uO)) , nO)) + E/-/ (nO-), BnO-)) 



Substituting n = n'"^ — X]f=i n(-'),we have 



E E c^P 



E E 



exp 



(z + u(0)) , (n - Eti' n(^)) ) + E/r^^ ((z + u^)) , n«) 



((z + u(")) , n) - Eti (("^"^ - ^^'^) : n^"^) 



n - Eti n« , S n - Eti n« ) + E'P (n(^), Sn^)) 



(43) 
(44) 



(45) 



(46) 



E E 

neZ3 n(i),,n(-'-i)« 



exp 



((z + u(0)) , n) - E;=/ ((u(°) - u(^) + Bn) , nW)] + 
n, Bn) + (Eti' n«, BEt/ n«\ + E/J/ (n^), Bn^) 



(47) 



^ exp ITT 2 (z + u(")^ , + (n, Bn) 



E ^^p 

n(i),.,n(J-i)eZ9 



Et'/((u^^'^-u^"^-Bn),nW) 



Et"/ n« , B E t'/ n« ) + E/=i (nW , Bn« ) 



(48) 



exp ITT 



2((z + u(°)) ,n) + (n,Bn) 



u(i) - u(0) - Bn 
u(j-i) _ u(o) _ 5n 



2B .. B 
2B .. 
B .. 2B 



(49) 
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The next substitution is n Jm + e, with n e , e e and thus X^nezs ~ Smezs EeiEZ" ' "^^^ 



standard properties of theta functions M,M that 



6l(z + Am|A) = cxp{-i7r[2(z,m) + (m,^m)]}6l(z|yl), if z eC", AeC"''" and m e Z" 
the expression becomes 



^ ^ exp|i7r 2<^ z + u(°) 

u(o) - u(i) + B ( Jm + e) 



, ( Jm + e)^ + (( Jm + e) , B ( Jm + e)> | 



u(0) -u^-'-i) +B(Jm + £) 



2B 
B 



B 

2B 



J2 ^eyi-p^^in ^2 (z + u(°)) + B ( Jm + e) , ( Jm + e)^ | 

u(o)-u(i)+Be 1 r 2B .. B 1 r m 
u(0) _ „(./-!) +Be _ B .. 2B _ _ m 

^ exp iTT /2 fz + u^"A + B ( Jm + e) , ( Jm + e] 



2B .. B 
B .. 2B 



X exp ( —ITT 



2 (E/=/ - u^^'' + ^e) . + (mJ (J - 1) Bm) 



u(0) _ u(i) + Be 
u(") - u(-'-i) + Be 



2B 
B 



B 

2B 



J2 exp {iTT [2 ((z + uo) , e) + (e, Be)]} (-^z+E/^o u^'' + -5e| JB 
u(o) _ u(i) + Be 



X 9 



u(o) _ u('^-i) + Be 



2B 
B 



B 

2B 



i.e. we have obtained the right hand-side of ( p^ ) under condition that Ej=o^'^^"'^ = 0- 



(50) 



(51) 



(52) 



(53) 
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7 APPENDIX 2 



Algorithm for derivatives of W-function. 

The VK-function derivatives and thus the combination of derivatives of r-function leading to some hier- 
archy of pde in N+l space-time can be derived algorithmically using e.g. the Mathcmatica program. 
Starting from the more useful form of ( p3| ) 

expg [inr (z + u^) - Inr (z)] ='£w (w^^), w(^-i); e) (54) 

with u^*^) = u^*^) (w^^\ w*^"'"^)) and differentiating with respect to the components of w^'^) vectors we 
obtain the expressions for combinations of derivatives of r -function logarithms. A simple algorithm for this 
calculation we illustrate here giving two examples. 



J ==3; j = exp(i27r/J) =exp(i27r/3); w^^^ e , fc=l,2 
Qp+q 

f)P+q 



W 



,(2) 



(55) 



exp 



L (w(i) + w(2)) + L (jw(i) -I- /w(2)) -f L (/w(i) + jw(2)) - 3L (0) 



J = 5; j = exp(i27r/J) = exp(i27r/5); w^*^) e fc = l,...,4 
Qp+q 



Qp+q 



exp 



(4)\ = 



.(1) 



L (j4w(i) 



F w(2) - 



f w(3) - 



5L (0) 



-j3w(4)) 
-j2w(4) 



(56) 



where have used a shorthand notation L (w) := Inr (z + w) . 
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